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Abstract. We report on experiments to measure the temporal and spatial evolution of packing arrange- 
ments of anisotropic and weakly confined granular material, using high-resolution 7-ray adsorption. In 
these experiments, the particle configurations start from an initially disordered, low-packing-fraction state 
and under vertical solicitations evolve to a dense state. We find that the packing fraction evolution is slowed 
by the grain anisotropy but, as for spherically shaped grains, can be well fitted by a stretched exponential. 
For a given type of grains, the characteristic times of relaxation and of convection are found to be of the 
same order of magnitude. On the contrary compaction mechanisms in the media strongly depend on the 
grain anisotropy. 

PACS. 45.70.-n Granular systems. - 45.70.Cc Static sandpiles; granular compaction - 81.05 Porous ma- 
terials; granular materials 



1 Introduction 

Granular media are not thermal systems: the thermal en- 
ergy fcflT is not involved in the evolution because it is neg- 
ligible compared to the variation of the gravitational en- 
ergy during the motion of a grain. Nevertheless they share 
with glasses a great number of properties - such as off equi- 
librium dynamics, aging and hysteresis. . . It has even been 
suggested that the relaxation of a granular medium under 
weak mechanical perturbations, such as shaking, has a for- 
mal analogy with the slow dynamics of out-of-equilibrium 
thermal systems [112] . This analogy is based on the as- 
sumption that the most important parameter in the sys- 
tem is the geometry and not the interaction between par- 
ticles or the driving energy. Consequently, although me- 
chanical agitation is neither stochastic nor isotropic, re- 
laxation of granular media is often presented as an ideal 
system to study out-of-equilibrium dynamics 0. 

The first experiments of slow relaxation in granular 
media submitted to external mechanical excitations have 
been carried out in Chicago | lloj . By implementing, at 
regular intervals, identical vertical taps of controlled in- 
tensity, they show that the packing fraction of the media 
increases slowly with the number of taps. This has been 
done for both spherical 0] and anisotropic grains [5]- For 
isotropic grains, the evolution of the strongly confined me- 
dia is correctly described by the inverse of the logarithm 
of the number of taps described in equation 

* = *~-(*--*° \ + BHi + t/Ty (1) 

The four parameters adjusted during the fit are <?o the 
initial packing fraction, <?oo the final packing fraction (for 



t — > 00), t, a characteristic time of the evolution and B a 
parameter without any physical interpretation. 

Other experiments have been carried out in quite dif- 
ferent systems: compaction under shearing |bl7| or under 
tapping but with a larger vessel 8 9 1011 1| - and thus 
a weaker confinement - than the one used in 0]. In [H] 
the best fit found is not the one used by the Chicago 
group but the so called KWW's fit (Kohlrausch, Williams, 
Watts) |12ll3j . This law is a stretched exponential (Eq.EJ: 

<P = $ 00 ~($oc~$a)eM-(t/T) P ), (2) 

where r is a characteristic time and (3 measures the slow- 
ing down compared to the simple exponential law. Here 
and <?oo are not fit parameters but experimental data: <Pq 
is the first point of the curve (the initial packing fraction) 
and the mean value of the packing fraction at the sta- 
tionary state. The tapping intensity is found to rule the 
characteristic time of relaxation according to an Arrhenius 
behavior relation. All these characteristics are similar to 
those obtained in strong glassy systems and thus these 
results confirm the analogy between glasses and granular 
media explained in the introduction. 

The aim of this article is to check the validity of the 
results obtained in |8I9| with non-spherical shaped grains. 
Indeed, to reduce the number of parameters, most pre- 
vious studies deal with spherical grains |4lbl7l8ll(JI9lll| . 
Nevertheless actual granular material are far from being 
perfectly isotropic. Grain anisotropy modify geometrical 
frustration and may change packing behavior and convec- 
tion during compaction. 

This paper is organized as follows. In the next sec- 
tion, the experimental setup is described. Results on the 
relaxation and on the dynamics of compaction are shown 
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in section 3. The last section contains a summary of the 
findings and conclusions. 

2 Experimental setup 

The experimental setup consists of a glass cylinder of di- 
ameter D w 10 cm filled with about 600 g of grains (corre- 
sponding to a height of roughly 10 cm). The method used 
to build our initial packing, already used for spheres [8I9| . 
is reproducible (we obtain the same initial packing frac- 
tion and the same packing fraction evolution for a given 
tapping amplitude). It consists of using two grids. The 
first one is placed at the bottom and inside the cylinder. 
The grains are then poured through the second one which 
is placed at the top the container. Once the filling proce- 
dure is finished, the first grid is pulled through the medium 
which is then dilated. This allows to obtain a low initial 
packing fraction. 

The container is placed on the plate connected to an elec- 
tromagnetic exciter (LDS V406) which induces a vertical 
displacement of the plate. The container is in this way 
shaken at regular intervals [At = Is) by vertical taps. 
Each tap is created by selecting an entire period of sine 
wave at a constant frequency / = 30 Hz. The resulting mo- 
tion of the whole system, monitored by an accelerometer 
at the bottom of the container, is however more compli- 
cated than a simple sine wave. At first, the system under- 
goes a positive acceleration followed by a negative peak 
with a minimum equal to —-f m ax as showed in Figure ^ 
Moreover, it can be seen on this figure that the acceler- 
ation created by the fall of the media on the bottom of 
the cylinder, —jfaii, may be different from —j ma x- Indeed 
we only used taps for which the media took off from the 
bottom of the glass cylinder. It means the tap intensity is 
always chosen upon the lift-off threshold. After the applied 
voltage stops, the system relaxes to its normal repose po- 
sition. Sequences of 10 4 to 10 6 taps are carried out with 
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Fig. 1. Example of dynamical signal monitored by the ac- 
celerometer for r ~ 3.5. The acceleration created by the fall of 
the media on the bottom —"/fail is different from the minimal 
acceleration applied to the system "/ max . 



an intensity r defined as follows: r — J max /g where g 
is the gravitational acceleration. Here we call "time" the 



number of taps and the "dynamics" is the succession of 
static equilibrium induced by the taps. 
The average volume fraction in the bulk <P is measured 
using the transmission ratio of 7-ray beam through the 
media. This setup allows to determine the vertical density 
profile too. The link between •P and the ratio N/Nq, where 
N and No are the activities counted respectively by the 
detector with and without media, is: <P w -hi(N/No). In 
this equation, \i is the absorption coefficient of the beads 
determinated experimentally |8l9j . The collimated 7-beam 
is nearly cylindrical with a diameter of 10 mm and inter- 
cepts perpendicularly the vertical axis of the cylindrical 
container. An acquisition-time of 60 s was found to be 
a good compromise between the intrinsic uncertainty of a 
radioactive beam and the total duration of an experiment. 
We then achieve a precision A$> — ±0.003. With the aim 
of limiting the duration of the experiments (from 20 hours 
to more than a month according to the total number of 
taps) and of avoiding redundant information due to the 
very slow evolution of the system, the measurements are 
spaced out in time, on a logarithmic scale, with 2 mea- 
sures of packing fraction profile and 50 measures of mean 
packing fraction of the sample ($) per decade (except 10 
for the first decade). 

The anisotropic grains used are rices of different shapes 
(see Fig. |2J): long grains (basmati rice) and short grains 
(round rice). Results obtained with spheres in previous 
works |8I9| are also reported. In order to quantify the grain 




Fig. 2. Picture of the grains used: basmati rice (left), round 
rice (center) and glass beads (right) 



anisotropy each grain is assimilated to an ellipsoid and the 
mean ratio between the two axis is measured. We obtain 
2.5 for the basmati rice and 1.5 for the round rice. It is 
worth noting that the size distribution is larger for the 
round rice (size dispersion around 20%) than for the bas- 
mati rice (less than 10%). 
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Fig. 3. Evolution of the packing fraction of a) basmati rice with F = 2.4, b) round rice with F — 2.4, c) basmati rice with 
r — 6 and d) round rice with F = 6. Chicago group's fit (dotted line) and KWW's fit (dashed line) are reported for each curve. 
For all the fits we used the Levenberg-Marquardt algorithm and the parameters #o an d $oo are set to the experimental values. 



3 Dynamics of compaction 

3.1 influence of the shape on the initial and final 
packing fraction 

For disordered packings of spheres, computer simulations 
as well as experiments have shown that the maximal pack- 
ing fraction is set by the so-called random close packing 
(RCP) w 0.64). Moreover the minimal packing frac- 
tion for a mechanically stable sphere packing (the random 
loose packing) is often set to <P m 0.55. Such estimations 
for disordered packings do not exist for anisotropic-grain 
packings. As mentioned above, the method used to build 
our initial packing is reproducible. It gives <P ks 0.56 for 
round rice and <P « 0.55 for basmati rice, which means 
that, in our experiments, the more anisotropic the grains 
are, the less important the initial packing fraction is. This 
can be explained by the fact that arches - and thus large 
pores - form more easily with anisotropic grains. Figure |3 
reports the evolution of the packing fraction with the num- 
ber of taps for the two kinds of rice and for r = 6 and 2.4. 
The first observations is that, whereas final packing frac- 
tions for round rice are comparable to what is obtained for 
compaction of sphere packings, they can be much higher 
for the basmati rice. Nevertheless we do not observe, con- 
trary to Villarruel et al , an ordered phase in which the 
grains align vertically . This order seems to be due to the 
strong lateral confinement. Indeed in 5 the length of the 



rods used is the same as the diameter of the cylinder so the 
steric constraints are huge. In our experiments the smaller 
ratio between the diameter of the cylinder and the length 
of the grain is about 25. Nevertheless, as explained below, 
we can observe for given values of r, an order which is 
different from that observed in [5]. 

3.2 Dynamics of compaction for basmati rice 

Depending on the value of r, two types of phenomena 
can be observed in basmati rice: for high values of tap- 
ping intensity r, convection is observed in the whole me- 
dia whereas for low values of r convection is localized in a 
part of the packing and is not steady. As it will be shown, 
compaction mechanisms strongly depend on the convec- 
tion and thus on r. It should be pointed out that we only 
consider values of r larger than the lift-off threshold jTH] . 
For high tapping intensity (typically r > 3) convection 
takes place within the whole medium. After about ten 
taps, two convection rolls appear but this situation be- 
comes unstable. One of the rolls progressively disappears 
and after this transient the whole medium is fulfilled by 
only one roll and the free surface of the medium is tilted 
from the horizontal (for example about 20° for r = 6). 
This can be observed at the vessel wall as showed in Fig- 
ure fj] With regular acquisition of images after each me- 
chanical perturbation, a characteristic time for the evolu- 
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Fig. 4. sketch of the convection obtained for basmati rice and 
for r — 6 (left) and snapshot of the media during relaxation 
(right) where convection roll can be clearly seen on the side 
walls 




Fig. 5. sketch of the convection obtained for basmati rice and 
for F = 2.4 (left) and snapshot of the media during relaxation 
(right) where the compact and ordered part can be clearly seen 
on the side walls 



tion of convection rolls, T conv , can be extracted. Indeed, 
we can track grain movements at the vessel sides and mea- 
sure how long it takes a grain to revolve around the center 
of the convection roll. For basmati rice submitted to taps 
of intensity r — 6 we find r conv w 50. The evolution 
of packing fraction during this experiment, is reported in 
Figure Ot- Chicago group's fit and KWW's fit (using the 
Levenberg-Marquardt algorithm) both correctly describe 
the experimental evolution and they give extremely simi- 
lar correlation coefficients (« 0.991). However, the evolu- 
tion of the former fit near the stationary state is slower 
than the experimental result. The fit provides the follow- 
ing values: t = 1150 and B — 165 for Chicago group's law 
and t = 30 and (3 — 0.58 for the KWW's law. The value of 
tkww is of the same order of magnitude than T conv and 
this remains valid for all our experiments done for r > 3. 
Note that for the fits the parameters <Poo and <Po are set 
to the experimental values. 

The behavior of our system is totally different for taps 
of low intensity (_T < 3). Two convection rolls roughly 
equally-sized and localized near the free surface appear. 
The other part of this packing (below the rolls) is not 
submitted to convection. As for the high tapping intensity 
case, this situation is unstable. Contrary to the previous 
case the two rolls do not merge: their size continuously 
decreases until they disappear. For basmati rice and for a 
solicitation of intensity r — 2.4, the convection is stopped 
in most of the packing after r ~ 1500. Convection rolls 
do not exist anymore, except near the air-media interface 
(see Fig.UJ). Another interesting point is that an ordering 



process is observed at the bottom of the vessel: the convec- 
tion rolls compact the lowest part of the packing and align 
the rice grains horizontally. This explains why the packing 
fractions obtained can be so high. This order creation is 
different from that observed by Villarruel et al [5] where 
the grains are preferentially oriented vertically. Thus the 
order observed in is created by the side walls whereas 
the order observed in the present work is created by con- 
vection. Note that the aspect ratio of the grains is proba- 
bly also an important parameter. The evolution of packing 
fraction during this experiment is showed Figure On 
this figure, it can be seen that the fit from Chicago group's 
law (logarithmic law) does not describe the experimental 
evolution very well. Indeed, Chicago group's law increases 
too fast before the steady state and reaches it too late. 
The parameters of Chicago group's law are r = 375 and 
B = 2 (with a correlation coefficient of 0.972). We re- 
call that for the fit we set the values of equal to the 
experimental final packing fraction. A better correlation 
coefficient can be obtained (0.996) if <Poo is a free param- 
eter but its value (0.86) is then unrealistic and well above 
the value of the plateau observed in Figure^. So for this 
law, the parameter r for r = 2.4 is smaller than r for 
r = 6 whereas the evolution slows down with the de- 
crease of r. This characteristic time which comes from 
Chicago group's fit does not show a physical behavior and 
does not describe our results. On the contrary KWW's law 
correctly fits the experimental evolution with the following 
parameters: r = 1380 and (3 = 0.31 (with a correlation co- 
efficient of 0.998). Like for the high- intensity solicitations, 
tkww the time extracted from KWW's fit and the char- 
acteristic time of the convection rolls are similar. So in all 
our experiment, tkww can be correlated with a real time 
of evolution. 

It should be pointed out that such kind of correlation is 
not possible in the experiments of Chicago because the 
use of a very narrow vessel prevents convection. Therefore 
the convection mechanisms are not the same in the two 
setup. 

3.3 dynamics of compaction for round rice 

For round rice, we did not observe the same change of be- 
havior for r ss 3. The convection always took place in the 
media, as for spherical beads |8I10| . And, as for spherical 
beads, KWW's fit describes better the experimental evo- 
lution of the packing fraction than Chicago group's law 
(see Figures |3Jd and ^) . However we could see a strong 
decrease of the convection speed for intensity taps under 
r « 3. More results about convection can be found in |14j . 

3.4 Dependence of the relaxation on r 

For our experimental setup, KWW's law is always the 
best fit. For all the grains used the parameter t follows an 
Arrhenius behavior (see Fig. Ep) 

r(r)= r expf^V (3) 
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The values of the parameters are ro « 1.06 and 2"b ~ 7 
for glass beads, 1.14 and 12 for short rice and 1.56 and 16 
for long rice. Even if the interpretation of To is not easy, 
the value of Zo is linked to the height of energy step in 
the phases space. Indeed, more anisotropic are the beads, 
larger is the distribution of energy in phases space and so 
greater must be the typical intensity of a tap to explore 
this space. 




Fig. 6. (a) Characteristic time of the stretched exponential 
as a function of the inverse of the tapping intensity r for glass 
beads (filled circles), round rice (open triangles) and basmati 
rice (filled squares) . (b) (3 exponent of the stretched exponential 
fits as function of the tapping intensity r for glass beads (filled 
circles) , round rice (triangles) and basmati rice (filled squares) . 
Glass bead data are taken from [H] and UJ. 




Fig. 7. Final packing fraction as a function of the tapping 
intensity F for glass beads (filled circles), round rice (open 
triangles) and basmati rice (filled squares). Glass bead data 
are taken from [Hj and UJ. 



4 Conclusion 

We have studied the effect of grain anisotropy on gran- 
ular compaction under vertical tapping, and weak lateral 
confinement. Using a 7-ray adsorption apparatus, we have 
reported the evolution of the packing fraction as a func- 
tion of the number of taps. We observe that the main 
features of granular compaction do not qualitatively de- 
pend on grain shape. In particular, the packing fraction 
evolution can be well fitted by a stretched exponential in 
all the cases. The fit parameters as well as the final pack- 
ing fraction have the same qualitative behavior and the 
fit parameter r is found to be correlated to the experi- 
mental characteristic time of the convection. Nevertheless 
we observed that the convection in the granular media, 
and thus the compaction mechanisms, strongly depend on 
grain anisotropy. 



The evolution of /3 is reported in FigureHJ). It is qual- 
itatively the same for the different kinds of grains. For a 
given value of tapping intensity Z 1 , the difference to a per- 
fect exponential (/3 = 1) is more important for anisotropic 
grains. The stretching of the exponential, and thus the 
coefficient is /?, is due to a wide characteristic time dis- 
persion. Since the range of metastable states is wider for 
anisotropic particles the characteristic times are more dis- 
persed and thus the coefficient (3 lower. 

The last quantity reported is the final packing frac- 
tion (Fig. EJ) . Although the general shape of the curves is 
similar, the quantitative measures depend on grain shape. 
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